Counting statistics of phase slips in superconducting nanocircuits 
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We perform measurements of switching current on different types of superconducting weak links 
and systematically study higher moments of the switching current distributions. We employ two 
types of devices in which a weak link is formed either by a superconducting nanowire or by a graphene 
flake subject to proximity effect. We find that, independently on the nature of the weak link, the 
third moment (skewness) of the distribution is universally close to —1 both in thermal and quantum 
regimes, thus implying manifestly non-Gaussian statistics. The fourth moment (kurtosis) also takes 
a universal value close to 5. The discovered universality of skewness and kurtosis is confirmed by 
the analytical model. Our numerical analysis shows that introduction of extraneous noise into the 
system leads to significant deviations from the universal values. We suggest to use the universality 
of higher moments as a robust tool to check against undesirable effects on noise. 

PACS numbers: 74.78.Na,74.50.+r,74.25.Sv,74.40.-n 

Introduction.-The field of quantum noise has recently seen rapid development caused both by its growing significance 
in many areas of condensed matter physics as well as by a constant improvement in the capabilities of high precision 
measurements (T]. Perhaps the most intensively studied question to date is related to the statistics of the charge 
transport in mesoscopic conductors. In such systems probability distribution of current fluctuations, the so-called full 
counting statistics (FCS), was rigorously derived for various normal and superconducting circuits [2] and tested in 
the state of the art measurements of the third moment (skewness) of current fluctuations [3] . Given that charge is a 
quantum mechanical conjugate variable to the phase, it is of fundamental interest to study corresponding statistics of 
phase fluctuations. Superconducting nanowires and related proximity devices offer a natural platform for this purpose, 
which we explore in the present work. 

The macroscopic quantum tunneling (MQT) of the phase across a current-biased Josephson junction [4 or a 
superconducting nanowire [5] is arguably the most profound and well-known manifestation of quantum fluctuations 
at the macroscopic level. This phenomenon is observed by registering phase slip events [6], which proliferate at currents 
close to the critical and drive transitions between supercurrent-carrying and dissipative parts of the I-V curve [THlOj. 
MQT is usually described in the framework of Caldeira-Leggett model [TT] that treats quantum or thermally activated 
transitions of the phase between neighboring minima of a tilted washboard potential in the presence of a dissipative 
environment (121 113] . Complimentary approaches employ an effective action for BCS superconductors [141 115] . 

Unlike experiments associated with the charge transfer where measurements of each moment of the FCS is beyond 
current experimental capabilities, experiments on switching current allow one to reconstruct full distribution of su- 
perconducting phase fluctuations since a single phase slip is sufficient to drive the system into resistive state [51 [TB] by 
creating a hot spot [171 118j . Thus there exists a one-to-one correspondence between phase fluctuations - phase slips - 
and switching current events. 

Here we report a systematic study of higher moments of the switching current distribution as a function of tem- 
perature and other parameters of our devices. The higher moments under investigation include skewness S that 
quantifies an asymmetry of the distribution, and kurtosis K that is a measure of its peakedness (for definitions see 
below). We present evidence, both experimental and theoretical, that these higher moments are in fact universal 
constants: S ~ — 1 and K sa 5. Surprisingly, the observed crossover from a classical escape mechanism (i.e., the 
thermal activation) to a quantum one (i.e., quantum tunneling from a metastable energy minimum) does not lead to 
any noticeable changes in these moments. We evince this universality using two types of samples, namely graphene 
junctions under the proximity effect as well as ultra-thin superconducting nanowires. Apparent universality of S 
and K has to be contrasted with the behavior of the switching current dispersion (second moment a) that exhibits 
nontrivial temperature dependence: the power-law [19j . a cx T 2 / 3 , in the thermal regime and a cx const in the quan- 
tum regime [SHTUl ITS] . We also demonstrate, numerically, that presence of extraneous noise leads to a substantial 
reduction of the universal moments. This observation provides an independent tool to assess the relevance of noise to 
interpretation of experimental data. 

-Devices.-Nanowire samples were prepared [101 I20j by depositing carbon nanotubes across a 100 nm wide trench 
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on a silicon chip, coated by a film of SiC>2 and a film of SiN. A film of 10-20nm of Mo7gGe24 was sputtered onto 
the chip, covering the top SiN surface and the nanotubes crossing the trench. Thus the suspended segments of 
nanotubes were converted into nanowires. Uniform wires were selected using SEM, and the MoGe film was patterned 
by photolithography, to define contact pads (electrodes). In such devices the selected nanowire serves as the only 
conducting link connecting the superconducting thin-film electrodes, positioned on the opposite sides of the trench. 
Importantly, there is no additional contact resistance between the nanowire and the contact pad since the wire 
transforms seamlessly into the pad while both are made in the same sputtering run. 

Graphene flakes were deposited onto an Si02 surface by mechanical exfoliation [21]. E-beam lithography was 
utilized to pattern the electrodes into a comb shape. After the resist was exposed and developed, we deposit, by 
means of thermal evaporation, a 4 nm Pd film (so-called sticking layer) and a 100 nm Pb film on the top. Lift-off was 
performed by placing the sample in an acetone bath for five minutes, sonicating it for ten seconds every other minute. 
The lOOnm Pb layer induces superconductivity in the graphene through the proximity effect. The samples were 
measured in a He-3 cryostat. Electromagnetic noise was filtered from the system using 7r— filters at room temperature 
and a copper powder and silver-paste radio-frequency noise filters at low temperatures. 

A sinusoidal bias current, having an amplitude greater than the critical current of the device, was applied across 
each sample. As the current increased from zero to its maximum, the voltage across the sample demonstrated a 
sudden jump from zero to some large, non-zero value, indicating the system switched from a superconducting state 
to a normal, resistive state. The value of bias current at which the jump took place was recorded as the switching 
current. Then the bias current returned to zero, and the system once again became superconducting. This process was 
repeated N = 10 4 times at each temperature for nanowire (NW) samples, and N — 5000 times at each temperature 
for superconductor-graphene-superconductor (SGS) samples. Each measurement gave slightly different value of the 
switching current, due to inherent stochasticity of the phase slips, thus producing switching current distributions. 
The skewness and kurtosis of each distribution was calculated from the recorded data by using standard expressions: 

N N 

S = ^2 {I s w,i~ (Isw)) 3 /p 3 and K = ^2 (I s w,i~ (Isw)) 4 /& A , where each I SWt i represents an applied bias current at which 

i=i ' i=i 
a switching event took place, (I sw ) is the mean switching current, and a is the standard deviation of the switching 
current distribution. 

Experimental results.-We first discuss the effect of tem- 
perature on the skewness and kurtosis of the switching cur- 
rent distributions. We find that in both superconducting 
nanowire samples (Figs. |Tp-|T]l) , and graphene proximity 
junctions (Figs. [2j:-[2jl) the skewness and kurtosis are con- 
stant with temperature. Surprisingly, these moments are 
identical within experimental uncertainty for the two quali- 
tatively different systems. The value of the skewness in both 
cases is near —1, and the value of kurtosis is near 5. In 
nanowire samples, these moments remain constant even as 
the system experiences a crossover from the high tempera- 
ture regime, at which phase slips are predominantly caused 
by thermally activation, to low temperatures, at which quan- 
tum tunneling of phase slips is responsible for the premature 
switching. This crossover is evident in Fig. [iji as the stan- 
dard deviation changes from the power-law at high tempera- 
tures to a constant value at low temperatures. The classical- 
to-quantum crossover temperature is typically in the range 
0.6-0.8 K for the studied samples (Fig. [I]). 

In SGS samples in addition to the temperature dependence 
we also study the effect of gate voltage V g on the skewness 
and kurtosis of the switching distributions. Both moments 
remain constant within the experimental uncertainty over a 
wide range of T and V g (see Figs. ^ and |Jl). It should 
be noted that, unlike NW samples, SGS junctions do not show crossover to the quantum tunneling dominated 
regime within experimentally tested temperatures. However, we do expect that such crossover might occur at lower 
temperatures, as recently reported [2"5|. 

Although most samples produced switching distributions with S ~ — 1 and K s=s 5, nanowire samples B and F have 
slight, but noticeable deviations from these values. Such deviations pull the moments closer to those of a Gaussian 
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FIG. 1: (a) standard deviation, (b) mean switching cur- 
rent, (c) skewness, and (d) kurtosis of the switching cur- 
rent distributions in nanowire samples A, B, C, D, E and 
F vs temperature T. The experimental values are repre- 
sented by the symbols. Simulation curves are shown by 
the solid lines. One point in figure (c) and two points in 
figure (d) lie outside the ranges shown. Fitting parame- 
ters used in simulation are summarized in Table I of of 
Supplementary materials |22| . 
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distribution. We attribute such trend to white noise present within the system which could be caused by e.g. depairing 
current fluctuations brought about by time-dependent morphological changes of our amorphous wires. Utilizing the 
same numerical model as used to fit the curves in Figs, [l] and [2] we examined how random noise in the critical current 
could affect the skewness and kurtosis of our nanowire samples. To do this, we allowed the critical current to vary 
about a mean value with Gaussian probability between the measurement runs. Results are shown in Fig. [3j We find 
that upon inclusion of such extraneous noise both skewness and kurtosis decrease in absolute values. Based on this 
observation we conjecture that reduction of the magnitude of higher moments from their universal values should be 
taken with caution and may signal about the presence of undesirable noise effects. 



Numerical simulations and fitting - The fitting curves pre- 
sented in Figs. [T] and [2] were built using the Arrhenius- 
type activation formula for the rate of phase slips (hereafter 

n = k B = i) MM 
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r(J,T) = fi(J,T) 



-U{I,T)/T 



-U(I,T)/T q 



(1) 



which accounts for both thermal and quantum escape pro- 
cesses. Here f2 is the attempt frequency, U is the energy 
barrier for the phase slip, T is the base temperature and T q 
is the quantum temperature used to model MQT regime ob- 
served in NW samples at low temperatures [25] • For both 
systems activation energy has power-law functional depen- 
dence on applied bias current 
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For SGS devices we took k — y/8/3 and rj 
and used critical current in the form 



(2) 

3/2 



Ic(T) = 



64ttT 



eRr. 



y 

" K + W n + y/2{W*+u n W n )X< 



A 2 (L/L n )exp{~L/L n ) 



(3) 



FIG. 2: (a) standard deviation, (b) mean switching cur- 
rent, (c) skewness, and (d) kurtosis of the switching cur- 
rent distributions in SGS sample Ills vs temperature at 
a gate voltages of IV, 3V, 5V and -IV. We use the same 
convention as in the previous figure. One point in figure 
(c) lies outside the range shown. Fitting parameters used 
in simulation are summarized in Table II of Supplemen- 
tary materials [22] . 



where Rn is the normal state resistance of a junction, A is 

the superconducting gap in the leads, oj n = (2n + 1)ttT, W n — \J A 2 + w 2 , L n = y / D/2uj n . The sum over n was 
taken until convergence (roughly 10 terms). Expression (|3| follows from the theory of disordered superconductor- 
normal metal-superconductor junctions [551 US]- It has to be stressed that ballistic theory of the proximity effect 
in SGS junctions [3U] fails to account for the temperature and gate voltage dependencies of the critical current for 
our devices (see Fig. pp). This observation is also consistent with the previous reports on the proximity effect in 
SGS systems [23l [26 | l3Tl - f36] . From the normal state resistance of our samples, we deduce typical mean free path 
I ~ 20nm, which correspond to the diffusion coefficient D ~ 50cm 2 /s. Because mean free path and Thouless energy 
Exh = D/L 2 ~ 80 jtteV are much smaller than the junction spacing of L ~ 300nm and energy gap A ~ lmeV, 
respectively, our SGS junctions correspond to long diffusive junction limit. 

For superconducting NWs there are two known models for U in Eq. If a wire forms a phase slip junction (PS J) 
then k = y/6/2 and v = 5/4 [HJ [37H39]. The corresponding expressions for k and r\ for the more thoroughly studied 
case of a Josephson junction (JJ) have the same values as above for SGS devices. It is worth noting that qualitatively 
the two models are very similar. Following the previous work |40) we model the critical current of NW devices by the 
phenomenological Bardeen's formula [H] 



/ c (T)=/ c (0)(l-r 2 /T c 2 ) 3 / 2 . 
Finally, for both SGS and NW systems the escape attempt frequency in Eq. ([T]) was described by 

n(i,T) = n (T)(i-i/i c (T)r 



(4) 



(5) 



with v = 1/4 for J J model, and v — 5/8 for PS J model. 

Eqs. ([l])-([5| were combined to determine the rate of phase slips. For a given set of parameters this rate was used 
to predict the switching distribution as a function of bias current and then calculate its mean, standard deviation, 
skewness and kurtosis. Such procedure was repeated at different temperatures to produce the temperature dependence 
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of the moments. In the case of SGS samples the above scheme was also repeated at different gate voltages. Parameters 
(f2o, / c (0), T c , T q ) for NW samples and {VLq, Rjy, T c and D) for SGS samples were then adjusted within the expected 
range of values until the predicted switching current and standard deviation vs temperature curves matched the data. 
These, along with the resulting skewness and kurtosis curves were used as fits to the data and are plotted as solid 
lines in Figs, [I] and [2] [221 H2] 

Analytical model-In this section we compute skewness and 
kurtosis by using an approach developed for the problem of 
escape from a metastable potential well subject to a steadily 
increasing bias field [TH1 H3] . We consider a general situation 
in which the phase slip rate of the system - either thermal 
or quantum - can be written in terms of reduced current i 
variable e = 1 — I/I c as 
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This form is general enough to cover all range of parameters 
relevant for our experiment on both types of devices. The 
powers a and b depend on whether the escape is quantum or 
thermally activated, while parameters A and B depend on 
the degree and type of damping (in particular, we estimate 
that our SGS junctions are moderately underdamped with 
the quality factor Q ~ 4). The distribution function for 
phase slips can be expressed in terms of the rate as 

P(e) = i-r(e)exp 



FIG. 3: Variation of the skewness (a) and kurtosis (b) as 
an auxiliary noise of the critical current is systematically 
added to the system. The horizontal axis in these graphs 
is the ratio of the standard deviation of the critical current 
of the wire and the standard deviation of the switching 
current which would occur if the critical current would be 
exactly constant in time. All these results are obtained 
by numerical modeling. 
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where |e| is a constant ramp speed. We are interested in central moments m„ of variable e i.e. moments defined around 
its mean value e: 



m n = ((e-e) n ) 



POO 

/ de(e-e)"P(e) 
Jo 



(8) 



where e = / °° deeP(e). Dispersion, skewness and kurtosis can be expressed in terms of central moments. To this 
end, it is convenient to introduce dimensionless parameter 



Z = In 



A/\e\ 



(9) 



which only weakly depends on the characteristics of the system in question. For self-consistency of the description 
this parameter should be large which can be achieved by tuning the ramp speed |e|. 

It is straightforward to show that moments of distribution Q can be written as an asymptotic power series in 
l/Z < 1 as follows: 



Q—nJb^njb 



Y,Z-ifi{nMZ) 



(10) 



Definition of the expansion coefficients fj are relegated to the Supplementary materials [53] because of their cumber- 
some form. Within the model fj depend on power exponents a and 6, and very weakly (as a double logarithm), on 
temperature-dependent parameters A and B, and the ramp speed |e|. This implies that temperature scaling of both 
moments ( e" ) and central moments ((e — e) n ) is fully dominated by the temperature scaling of B which is proportional 
to the height of the phase slip barrier {(e™), ((e — e) n }} cx B~ n / b {T). To determine proportionality coefficients one 
needs to use the explicit form of fj. Up to the order l/Z first two moments are given by 



- e=z i/b B -i/b (l+ V J* 
a 2 = m 2 = Z 2 ' h - 2 B- 2 ' b 



(eS + z¥ [ a7rV3 + (1 " 6)( " V3 ~ vy ' (1)) ] ) 



(11) 
(12) 
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We have defined v = (a/b) \nZ + 7, where 7 ~ 0.577 is the Euler-Masheroni constant, and ~ —2.404 is the 

tetragamma function |44j . Despite increasing complexity of the calculation the leading term in the third and fourth 
central moments are given by simple expressions: 

m^B-^Z^f-^l+Ss), (13) 



b 3 

7714 



= + (14) 



The correction terms are of the order {^3, 04} oc Z 1 . For example 83 = 60 l zhi [90a7r 2 t> (v — 1) — ll7r 4 (6 — 1) — 
180V'"(l)(a — v(b — 1))] [22 . We are now in the position to compute skewness and kurtosis and thus find: 

S = ~m 3 /m 3 2 /2 = 6 3/2 V"(l)A 3 + OiZ- 1 ), (15) 
K = nu/ml = 27/5 + 0(^ 1 ), (16) 

which are central results of this section. Remarkably, to the leading order in 1/Z -C 1, both skewness and kurtosis 
are given by universal numbers S ~ —1.139 and K s» 5.4, which are independent of the parameters of the system and 
are the same for both thermal and quantum phase slips. The magnitude of the correction terms 03, 04 is analyzed for 
different values of the ramp speed and different models of a weak link in Appendix |22j . 

Conclusion.-We have experimentally demonstrated the universality of higher moments - skewness and kurtosis - 
of the switching current distribution in superconducting nanocircuits. Our results are supported both by analytical 
modeling and by numerical simulations. We have also pointed out that the universality of higher moments is affected 
by extraneous noise and suggested to use this observation to detect the presence of unwanted noise in the data. 

This material is based upon work supported by the U.S. Department of Energy, Division of Materials Sciences 
under Award No. DE-FG02-07ER46453, through the Frederick Seitz Materials Research Laboratory at the University 
of Illinois at Urbana-Champaign. V. V. was supported by the U.S. Department of Energy, Division of Materials 
Sciences and Engineering under Award DE-FG02-08ER46544 and by the Theoretical Interdisciplinary Physics and 
Astrophysics Center at JHU. A. L. acknowledges support from Michigan State University and in part ARO through 
contract W911NF-12-1-0235. 

APPENDIX A: EXPANSION COEFFICIENTS fj(n,\nZ) 

Introducing new variables r — n/b and v = (a/b) I11Z + 7, we find that expansion coefficients introduced in Eq. (10) 
of the main text are given by the following expressions: 

f 1 (n,\nZ)=vr, (17) 

f 2 (n, \nZ) = a b (v + l)r + (v 2 + , (18) 

or 



l r „ , , 7T \ ar(r — 1) ( 9 
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r(r — l)(r — 2) / 3 n 2 v 



h{nMZ) = ~ [v 2 + l[2v(b-a) + 2b-3a} + ^\ + - v fe ^ lv 2 +v + (> 
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^"(1) , (19) 



, / , rys ar ( 3 k 2 v 2 ,„._,. 3(26 - 3a) / 2 tt 2 \ 3av(3b - a) 3(26 - a) . 2 2 . 
U(n, lnZ) = - (v 3 + — - ^'(1 + [v 2 + — K — 2 + - V - 7 + I 2 ) 



36 V 2 r w 26 V 6 / b2 b 

r(i — 1) (a 2 ( 2 7r 2 \ a / 3 n 2 v „ 2(6— a) / 2 7r 2 \ v(2b — 3a) 



(« + l) 2 + -5- « 3 + V-<( 1 ) 



6 2 V 6 y 26 V 2 v w 6 V 6 
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APPENDIX B: ANALYSIS OF CORRECTION TERMS S 3 AND 5 4 

For our results to be valid the ramp rate |e| must be small (otherwise the system will not be able to maintain an 
equilibrium from the initial state). The typical range of the sweep rates used in the experiments is |e| ~ (10~ 7 — 
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FIG. 4: Correction to skewness as a function of exponent a for different values of the ramp speed encoded by expansion 
parameter Z. Z = 15,20,25,100 from the bottom to the top curve. Notice that for most models summarized in Table III 
parameter a falls in the range where 83 is nearly independent of a and is determined solely by the ramp speed. In particular, 
for PS J model a — 3/8 and for J J model a — —1/4. 



10 _5 )s _1 . At operating temperature T ~ 0.3K the typical value of the escape frequency translates into the parameter 
A ~ 10 _12 s _1 . For the critical current at that temperature we have I c ~ 1CP 5 A which translates into the parameter 
B ~ 10 3 that gives the scale for the height of the potential barrier for the escape. This together leads to an estimate 
of the expansion parameter Z ~ 15. 

Now we present estimates that give some measure of the importance of correction terms S3 and 64 in our expressions 
for skewness and kurtosis versus power a and expansion parameter Z for the JJ model with b = 3/2. Given that 
theoretical calculations of the attempt frequency il (the pre-exponential factor in the escape rate) are quite challenging, 
the power-law exponent v (or, in the reduced notation of Eq. 6 parameter a) is uncertain and is still a subject of 
debate. To the leading order the expression for £3 reads 

S3 = ^4 [ 90a7r M^ - 1) - ^(b - 1) - 18(h//'(l)(a - v(b - 1))]. (21) 

From Fig. 1 we conclude that the magnitude of S 3 is about 10% of the universal value of S. The expression for S4 is 
implicit in Eq. (4) for f 4(12, hi Z) and also f§(n,hiZ). Taking a — for simplicity we find 

64 = 5Z^ [37?r4 " 157r ^"W " 10 ^ (4) (1)]- (22) 

For the value of Z estimated above (Z ~ 15) and for the JJ model with b = 3/2 one finds 64 « —0.679 which changes 
K by about 10%. These estimates provide a firm assurance for the consistency of our results. 



APPENDIX C: FITTING PARAMETERS 



TABLE I: Fitting Parameters for nanowire samples 



Sample Oo 




x ) J e (0) (M) T c (K) T q (K) 


A 


0.8 


11.02 


5.70 0.784 


B 


0.16 


12.01 


5.48 0.781 


C 


0.32 


13.27 


5.00 0.814 


D 


0.16 


9.38 


5.09 0.087 


E 


0.16 


4.25 


3.24 0.521 


F 


0.02 


5.46 


4.57 0.711 
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TABLE II: Fitting Parameters for SGS junction sample Ills 



Gate Voltage (V g ) flo (ns" 1 ) R N (Q) T c (K) D (cm'/s) 



1 


1.5 


138 


7.2 


51 


3 


1.5 


88 


7.2 


54 


5 


0.2 


72.5 


7.2 


57 


-1 


0.6 


210 


7.2 


48 



TABLE III: Parameters of the thermal and quantum escape rates for various degrees of Ohmic damping encoded by the quality 
factor Q of a junction. The table is compiled from previous studies of phase slip junction (PSJ) and Josephson junction (JJ) 
models [EHI5I E7H331 131 SSI US] • U stands for U(I = 0, T). 



Escape 


Damping 


A 


B 


a b 


Thermal-JJ 


Low 


Wn U /5irQT 


U /T 


1 3/2 


Thermal-JJ 


Moderate 




U /T 


-1/4 3/2 


Thermal-JJ 


High 


«0Q/27T 


U /T 


3/2 


Thermal-PCJ Moderate 


VUo/T 


Uo/T 


3/8 5/4 


Quantum 


None 


v /216t/ofi /7r 


36(7 /5«o 


5/8 5/4 


Quantum 


High 


V3£/ fio/Q 7 


3nU /QQ 


1 
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